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Abstract 

Ocean acidification (OA) has important implications for the persistence of coral reef ecosystems, due to potentially negative 
effects on biomineralization. Many coral reefs are dynamic with respect to carbonate chemistry, and experience fluctuations 
in pC0 2 that exceed OA projections for the near future. To understand the influence of dynamic pC0 2 on an important reef 
calcifier, we tested the response of the crustose coralline alga Porolithon onkodes to oscillating pC0 2 . Individuals were 
exposed to ambient (400 |iatm), high (660 |iatm), or variable pC0 2 (oscillating between 400/660 uatm) treatments for 14 
days. To explore the potential for coralline acclimatization, we collected individuals from low and high pC0 2 variability sites 
(upstream and downstream respectively) on a back reef characterized by unidirectional water flow in Moorea, French 
Polynesia. We quantified the effects of treatment on algal calcification by measuring the change in buoyant weight, and on 
algal metabolism by conducting sealed incubations to measure rates of photosynthesis and respiration. Net photosynthesis 
was higher in the ambient treatment than the variable treatment, regardless of habitat origin, and there was no effect on 
respiration or gross photosynthesis. Exposure to high pC0 2 decreased P. onkodes calcification by >70%, regardless of the 
original habitat. In the variable treatment, corallines from the high variability habitat calcified 42% more than corallines from 
the low variability habitat. The significance of the original habitat for the coralline calcification response to variable, high 
pC0 2 indicates that individuals existing in dynamic pC0 2 habitats may be acclimatized to OA within the scope of in situ 
variability. These results highlight the importance of accounting for natural pC0 2 variability in OA manipulations, and 
provide insight into the potential for plasticity in habitat and species-specific responses to changing ocean chemistry. 
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Introduction 

Anthropogenic C0 2 emissions have been rising steadily since 
the Industrial Revolution with the increase in fossil fuel 
consumption, agriculture, and land development [1], resulting in 
greater concentrations of C0 2 entering the surface ocean [2]. 
Ocean acidification (OA) refers to the equilibration of atmospheric 
C0 2 with the surface ocean, and the subsequent changes to 
seawater carbonate chemistry. These reactions are manifest as a 
decrease in mean ocean pH, carbonate ion concentration (CO3 ) 
and CaCO^ saturation state (£2), an increase in pC0 2 and 
bicarbonate ion concentration (HC0 3 _ ), and no change in total 
alkalinity (Ax). Saturation state is an important factor that 
influences rates of biomineralization, where CaCO :j precipitation 
is favored when Q is >1 [3]. Changes in Q. associated with ocean 
acidification have important implications for biogenic calcifiers 
and for ecosystems that are structurally dependent on carbonate 
platforms, such as coral reefs [3]. The changes in carbonate 
chemistry due to OA have raised widespread concern, particularly 
because they are occurring at a rate that is at least 2 orders of 
magnitude faster than during previous glacial periods [4]. The 
accelerated rate of OA may outpace the ability of marine 
organisms to adapt to the changing environment [5]. 



Near future projections for atmospheric CO z gas concentra- 
tions, such as in the Representative Concentration Pathways 
(RCPs) [6], are based on open ocean conditions where carbonate 
chemistry varies predictably over time and mean ocean pH is 
reported to be 8.02-8.07 [7]. Seawater chemistry in near shore 
coastal ecosystems, such as upwelling regions [8,9], coral reefs 
[10-12], estuarine systems, kelp forests [7] and natural C0 2 vents 
[13-16], is both spatially and temporally dynamic within a habitat 
[17]. Swings in diel pH frequently exceed the forecasted 
projections of ocean acidification for the end of the century 
[7,12,18-20]. The recent development of autonomous sensors that 
record high frequency in situ pH data over time [21,22] have shed 
light on the extent of pH variability within habitats and across 
ecosystems [7,12]. 

In coral reef ecosystems, particularly shallow reef flats, 
carbonate chemistry varies over the course of a day [17], with 
tidal cycles [23,24], and even between upstream and downstream 
habitats [25,26]. The shallow reef benthos often is dominated by 
biogenic calcifiers that can create high pC0 2 and low pH 
conditions by depleting total alkalinity (A x ) [27] and releasing 
C0 2 during the process of calcification [28,29] . The interplay of 
reef metabolism and carbonate chemistry thus creates a feedback 
that influences the magnitude of pH variability that resident 
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organisms experience [12,30-32]. A feedback also exists in the 
balance of photosynthesis and respiration, with peaks and lows in 
photosynthesis during the day and night respectively [10,33,34]. 
Exposure of resident marine organisms to dynamic swings in diel 
pH may influence the plasticity of organismal responses to 
changing carbonate chemistry, yet variation in pH has been 
incorporated in only a few ocean acidification manipulations 
[17,35-37]. 

To date, much of the ocean acidification literature has focused 
on the effects of decreased pH/increased pC0 2 on biological 
responses of marine organisms in static laboratory mesocosm 
experiments by maintaining steady C0 2 treatments over time 
[38]. Although the magnitude of responses varies, the overall effect 
of ocean acidification on marine calcifiers primarily has been 
negative [38-40]. These studies have been critical to advancing 
the understanding of C0 2 effects on marine organisms, but may 
underestimate the importance of natural fluctuations in pC0 2 on 
organismal responses. For instance, coral recruits exposed to 
ecologically relevant oscillations between ambient and high pC0 2 
experienced enhanced growth and survival when compared to 
recruits exposed to ambient or high pC0 2 conditions alone [36]. 
By understanding the role that natural variability plays towards 
influencing species responses to changing carbonate chemistry we 
can gain insight into the potential for organisms to acclimatize and 
adapt to ocean acidification. This issue is of primary concern 
because the accelerated rate of ocean acidification is unprece- 
dented over geological history [41] and may preclude the adaptive 
potential of some marine species. Little is known about species 
acclimatization potential to ocean acidification, and studies are 
beginning to explore the concept of acclimation in the context of 
ocean acidification [42]. Understanding how organisms already 
may be acclimatized in a rapidly changing environment can 
improve our understanding of how coral reef community structure 
may change in the near future. Marine ecosystems in which 
organisms are exposed frequendy and consistentiy to environ- 
mental fluctuations, such as shallow coral reef flats [7,10,12], 
provide an in situ platform in which to explore the potential for 
acclimatization to ocean acidification. 

Although much of reef net CaC0 3 deposition occurs through 
biogenic calcification by scleractinian corals [43], another 
important group of reef builders are crustose coralline algae 
(CCA) [44,45]. Crustose corallines serve many important ecolog- 
ical functions by contributing to reef primary productivity [46-48], 
facilitating settlement of coral larvae [49-51], and maintaining the 
structural integrity of the coral reef framework by cementing reef 
fragments together [52,53]. Corallines precipitate the most soluble 
polymorph of CaC0 3 , high-magnesium (Mg) calcite [54], and 
have been shown to be sensitive to ocean acidification. Exposure 
to high pC0 2 conditions reduces coralline calcification and growth 
[40,55-60], and inhibits the function of the chemical cue produced 
by corallines and/or their microbial assemblages that facilitates 
coral settiement [61-64]. Despite their important role in coral reef 
ecosystems relatively few studies have addressed the effects of OA 
on tropical coralline algae. Furthermore, studies have yet to 
explore the effects of variable pC0 2 on crustose corallines, or their 
potential for acclimitization. Porolithon onkodes (synonymous with 
Hydrolithon onkodes [65]) is an important reef-building coralline alga 
that is abundant throughout the Indo-Pacific in shallow reef 
systems and has been shown to promote settlement in coral larvae 
[49,51]. By exploring the response of P. onkodes from a naturally 
dynamic shallow reef flat to oscillating pH, we may better 
understand the adaptability of an important reef-building species 
to ocean acidification. 



The aims of this study were, 1) to determine the effect of 
oscillating pC0 2 on crustose coralline calcification, photosynthe- 
sis, and respiration, and 2) to determine the role of prior exposure 
to in situ pC0 2 variability in influencing coralline responses to 
ocean acidification. To accomplish the first goal, we conducted 
tightly-controlled laboratory mesocosm experiments in Moorea, 
French Polynesia, in which we exposed P. onkodes to ambient, high, 
and oscillating pC0 2 at levels ranging from ~400 u.atm to 
660 uatm, the natural range in pC0 2 on the back reef of Moorea. 
We hypothesized that exposure to high pC0 2 alone would 
decrease P. onkodes calcification, as shown in previous studies, and 
that exposure to pC0 2 oscillating within an ecologically relevant 
range would modulate the negative effect of high pC0 2 on 
calcification. Further, we hypothesized that exposure to high 
pC0 2 would enhance net photosynthesis in all corallines due to 
the potential promoting effects of increased dissolved C0 2 on 
carbon fixation. To explore the potential for acclimatization to 
ocean acidification conditions, we used P. onkodes individuals 
collected from either low pH variability sites (reef crest, upstream) 
or high pH variability sites (200 m shoreward of the reef crest, 
downstream) within the same reef habitat. W e hypothesized that 
individuals from downstream sites would calcify more in the 
variable pC0 2 treatment than those from upstream sites because 
they were acclimatized in situ to variable pC0 2 . To our knowledge, 
this study is the first to quantify the response of coralline algae to 
variable pC0 2 and to explore the potential for acclimatization of 
an important reef builder to high pC0 2 . 

Materials and Methods 

Site Description 

This study was conducted during May and June 2010 in 
Moorea, French Polynesia as part of the Moorea Coral Reef Long- 
Term Ecological Research site (MCR LTER) under a Protocole 
D'accueil (Scientific Research Permit) to RCC from the Delega- 
tion a la Recherche de la Polynesie Francais (unnumbered). 
Moorea is located 20 km west of Tahiti and is surrounded by a 
barrier reef that encloses a lagoon. Water circulation is driven by 
oceanic and wind-driven swells that break on the reef crest and 
force mostly unidirectional flow of water across the back reef and 
into the lagoon where it then exits through one of the reef passes 
[66]. In the back reef on the north shore of Moorea, water 
residence time is dependent on offshore wave height. The transit 
time for water flowing over the reef crest to habitats 200 m 
downstream is approximately 30-45 minutes [66] and is 
dependent on flow speed. The carbonate chemistry of water 
flowing across the back reef is altered by reef metabolism, with 
photosynthesis and calcification decreasing pC0 2 over the reef 
during the day, and respiration and dissolution increasing pC0 2 
over the reef during the night [29] . The farther a parcel of water 
travels from the source of oceanic inflow over the reef crest, the 
longer the residence time and exposure to reef metabolism and the 
greater potential for change in the carbonate chemistry. As a 
result, reef organisms inhabiting upstream locations (e.g., fore reef 
and reef crest) are exposed to relatively constant pC0 2 of 
incoming oceanic seawater, while the pC0 2 experienced by 
organisms living in downstream habitats varies on a diel basis and 
with the environmental conditions of light and water flow that 
modulate reef metabolic processes and determine water residence 
time over the reef [25,26,31,67,68]. 

Sample Collection 

Samples of the common crustose coralline alga, Porolithon 
onkodes, were collected from upstream habitats on the shallow fore 
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reef at a depth of 2-3 m on the north shore reef between Cook's 
and Opunohu Bays. Samples of the same species also were 
collected from the back reef approximately 200 m downstream of 
the reef crest at a similar depth. These two locations experience 
approximately the same flow [66], nutrient, and sedimentation 
[69] conditions. Individuals collected from the upstream location 
represent organisms exposed to a history of low pH variability, and 
individuals from the downstream location a history of higher pH 
variability [26,31]. Crustose coralline samples were collected using 
a diamond-grit hole saw (40 mm outer diameter, 36 mm inner 
diameter) attached to a pneumatic drill. Coralline algal thalli were 
separated from the underlying dead coral skeleton to yield 
approximately 0.5 cm thick cores of P. onkodes. Coralline algal 
samples were returned to the lab and carefully cleaned of 
epibionts. The bottoms of the algal disks were covered with 
marine epoxy (Z-spar) to prevent exposure of the underlying 
carbonate to potential dissolution in treatment conditions. Samples 
were kept in a water table with a high exchange of fresh seawater, 
under ambient light, for three days to allow for healing. Only 
healthy appearing samples without obvious tissue damage were 
used in the experiment. 

After samples were collected, the taxonomic identification of 
Porolithon onkodes was assigned based on Adey et al. (1982) [70] and 
Payri et al. (2000) [71]. We used gross morphological, reproduc- 
tive, and trichocyte characteristics to identify each sample using 
simple laboratory procedures. Specimens were dried at 60°C for 
48 hours, fragmented with diagonal cutters, and internal charac- 
teristics of the fragmented edges, including conceptacle size and 
shape, trichocyte arrangement, cellular organization, and margin 
morphology were examined under a dissecting scope. Identifica- 
tion of samples was verified by comparison to samples (N = 3) that 
had been previously confirmed by Dr. Robert Steneck (personal 
comm.). 

C0 2 Enrichment Treatments 

The C0 2 enrichment treatments were established in a 
mesocosm facility at the Richard B. Gump South Pacific Research 
Station consisting of 6, 150-L tanks, each with independent 
lighting (250 W metal halide fixtures) and temperature control. 
Mesocosm tanks were operated in an open configuration with each 
tank receiving fresh, unfiltered seawater pumped from Cook's Bay 
at a rate of 140 ml min -1 tank -1 resulting in approximately two 
complete turnovers of water each day. Each tank was fitted with a 
clear plexiglass lid to reduce atmospheric gas exchange and 
stabilize treatment pC0 2 levels. Light levels (PAR) were adjusted 
to maintain the light level in each tank at ~600 umol quanta m 2 
s 1 (on a 12:12 h lightdark cycle) that resulted in a daily 
integrated 12-h PAR similar to that measured at the collection 
sites (~26 mol quanta m~ 2 d -1 ). The temperature control for 
each tank was set to a target value of 28.5°C which represents the 
ambient temperature on the reef during the austral winter (May/ 
June). 

To simulate the effect of atmospheric CO2 on ocean surface 
waters, treatment tanks were bubbled with either ambient air or 
C0 2 -enriched air to establish targeted pC0 2 levels. C0 2 -enriched 
air was created using a solenoid-controlled gas regulation system 
(Model A352, Qubit Systems) that mixed pure C0 2 with ambient 
air to achieve the desired partial pressure (sensu Edmunds et al. 
(2011)) [72]. The targeted pC0 2 value for the high pC0 2 
treatment was 650 u,atm, which represents the pCO z at down- 
stream sites just before dawn and the projected atmospheric pC0 2 
in 50 years based on carbon emission scenarios [6]. The three 
pC0 2 treatment levels used in the experiment were ambient 
(—400 |J,atm), high pC0 2 (—660 |latm), and oscillating pCQ 2 



(samples were exposed to a treatment that alternated between 
400/660 natm). 

The physical parameters of each tank were monitored 
throughout the experiment including light, salinity, temperature, 
and carbonate chemistry. Temperature was measured daily in the 
morning, afternoon, and evening using a digital thermometer with 
an accuracy of 0.05°C (Traceable Digital Thermometer, Thermo 
Fisher Scientific). Light was measured daily beneath the surface of 
the water and at the sample locations (— 1 5 cm above the bottom 
of the tank) in each tank using a 2ti LiCOR quantum sensor 
connected to a LiCOR 1400 meter. Salinity was measured daily 
on water samples from each tank using a desktop YSI 3100 
conductivity meter, and pH (total scale) was measured daily using 
a spectrophotometric technique (described below). 

Seawater Chemistry 

Water samples were collected daily at 0900 hours in 20 mL 
glass scintillation vials for spectrophotometric analysis of pH (total 
scale). pH measurements were made immediately following 
collection using a spectrophotometer (Shimadzu UV-2450) and 
the indicator dye m-cresol purple according to standard operating 
procedure (SOP) 6 [73]. The temperature of the sample was 
controlled with a constant-temperature cell holder (Shimadzu) that 
maintained cell temperature at 25°C. The accuracy of the 
spectrophotometric technique was assessed by conducting pH 
measurements on certified reference material with a known pH 
(Tris Buffer in synthetic seawater, A. Dickson, Scripps Institution 
of Oceanography), and the mean accuracy of this technique was 
±0.09% (N = 16). Water samples were collected every two days at 
0900 hours for determinations of total alkalinity (A T ) by siphoning 
treatment water into 250-ml borosilicate sample bottles [73]. 
Water samples for A x were analyzed within 6 hours of collection 
after the water sample had equilibrated with laboratory room 
temperature. Total alkalinity (A T ) determinations were made using 
a Mettler T50 automated titrator to conduct modified, open cell 
potentiometric titrations [74] according to standard operating 
procedure (SOP) 3b [73]. Titrations were conducted on 100-ml 
samples at a temperature of 24°C (samples were equilibrated with 
room temperature) with a DG1 15-SC pH probe (Mettler Toledo) 
calibrated using a three-point calibration with NBS buffers 
(pH 4.00, 7.00, 10.00). The mean accuracy of the AT determi- 
nations was ±0.5% (N = 6) compared to Certified Reference 
Materials (Reference Material for Oceanic C02 Measurements, 
Batch 99, A. Dickson, Scripps Institution of Oceanography). 

Water samples were collected at each site over five days in order 
to characterize the carbonate chemistry at the upstream and 
downstream sites. Integrated water samples were collected during 
the day (from 0800 to 1800) and during the night (from 1900 to 
0400) by continuously sampling water (4 mL min ') with a 
submersible autonomous water sampler. Upstream and down- 
stream water samples were collected simultaneously. Water 
temperature was recorded at the time of sample collection and 
samples were returned to the lab and analyzed as described above 
within 24 h of collection. The results from all samples were 
processed in Microsoft Excel [74], and A T , pH, salinity and 
temperature were used to determine the remaining carbonate 
parameters using C0 2 SYS [75] with the dissociation constants of 
Mehrbach et al. (1973) [76] refit by Dickson & Millero (1987) [77]. 

The environmental conditions in the established treatments 
were stable over the 14-d experiment. Temperature, salinity, light, 
and total alkalinity were similar across all treatments and also 
similar to ambient levels in the field (Table 1). The pC0 2 
treatments created by bubbling ambient and C0 2 -enriched air 
resulted in distinct, non-overlapping pH and pC0 2 levels 
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Table 1. Carbonate chemistry measurements. 





Treatment 


Temp (*C| 


Salinity (psu) 


Light 1 


PH* 


TA (umol/kg) 


pC0 2 ((latml 


fie, 




Ambient 


28.7±0.02 


35.5 ±0.02 


607±15 


8.041 ±0.007 


2311 ±8 


417±16 


5.58 ±0.1 3 


3.72±0.09 


Ambient 


28.5 ±0.02 


35.5 ±0.03 


631 ±15 


8.035 ±0.005 


2312±7 


404±21 


5.66 ±0.1 5 


3.78±0.10 


Elevated 


28.5±0.02 


35.5 ±0.04 


651 ±17 


7.863±0.014 


2315±8 


667 ±30 


4.12±0.15 


2.75±0.10 


Elevated 


28.7±0.02 


35.5 ±0.03 


654±15 


7.863±0.013 


2314±3 


663 ±42 


4.17±0.19 


2.78±0.12 


Var/Amb 


28.7±0.02 


35.5 ±0.02 


668±16 


8.027 ±0.005 


2329±6 


432±6 


5.52 ±0.06 


3.68±0.04 


Var/Elev 


28.4±0.02 


35.5 ±0.02 


634±12 


7.865 ±0.01 4 


2322±3 


666±67 


4.17±0.32 


2.78±0.21 





Values (means ± SE) for physical variables in each tank over the course of the 14-d experiment. N = 14 sampling days for temperature, salinity, light, and pH, and N = 5 

sampling days for TA, pC0 2 , Q C a* an d H Ar . 

1 Light = photosynthetically active radiation (PAR, fimol quanta), 

"'"pH = total scale pH units, Q Ca = the saturation state of calcite, Q Ar = the saturation state of aragonite 
doi:1 0.1 371 /journal.pone.0087678.t001 



(Figure 1A) and saturation states of calcite (Qcalc) (Table 1). The 
pC0 2 levels in the elevated treatments (680±15, mean ± SE) 
were 65% higher than in the ambient treatments (411 ±14) 
(Figure 1A) and Q (:a i c was reduced by 28% from 5.67 ±0.04 to 
4.09 ±0.06 as a result of a reduction in the concentration of 
carbonate. There was no temporal variation in pC0 2 within the 
treatment tanks over the 2, 24-h periods that were sampled 
(Figure IB). 

Experimental Incubations and Response Variables 

Algal samples were allocated randomly to one of three pCO a 
treatments: 1) ambient pC0 2 (~400 (iatm), 2) elevated pC0 2 
(~660 (iatm), and 3) a variable pC0 2 treatment that cycled 
between 400 and 660 |j,atm (Figure 1). The elevated pC0 2 was 
chosen based on pre-dawn measurements of pC0 2 at the 
downstream collection site (Table 2). Prior to placing coralline 
algal samples into the treatments, each sample was buoyant 
weighed with the epoxy base (to the nearest mg) [78]. Six samples 
from each habitat origin (upstream, downstream) were placed on 
elevated PVC racks in each tank (15 cm above the bottom). The 
ambient controls and pC0 2 treatments were fully replicated twice 
(N = 2 tanks). The variable pCO z treatment was accomplished by 
placing samples in an ambient pC0 2 tank for 18 h (from 0600 to 
0000) and then transferring the samples to a dedicated elevated 
pC0 2 tank (i.e., samples in different treatments were not mingled) 
for 6 h (from 0000 to 0600) each day to simulate the timing of 
elevated pC0 2 in downstream habitats in the field (Figure IB). At 
each time of transfer, samples from the ambient and elevated 
treatments also were removed from their tanks and returned to the 
same tank to serve as a sample handling control. 

After 14 days in treatment conditions, all coralline algal samples 
were buoyant weighed again and 3 samples from each origin and 
pC0 2 treatment replicate were chosen randomly for measure- 
ments of photosynthesis and respiration. Photosynthesis measure- 
ments were estimated during the day from changes in dissolved 
oxygen during 30-min incubations of each sample in a 102-ml 
acrylic chamber fitted with a PreSens dipping oxygen optode 
connected to a Fibox 3 transmitter (Precision Sensing GmbH, 
Germany). The oxygen probe was calibrated at the start of 
incubations using a 2-point calibration with temperature correc- 
tion in water-saturated air (100%) and seawater with no oxygen 
(supersaturated sodium dithionite, Na 2 S 2 0 4 ). Oxygen concentra- 
tions automatically were corrected for temperature using a 
PreSens probe that simultaneously measured temperature in the 
incubation chamber (Precision Sensing GmbH, Germany). A 



water jacket surrounding the acrylic chamber was connected to a 
circulating water bath (Lauda) that provided temperature control 
(28°C). A stir bar at the base of the chamber created vigorous 
water motion. Light was provided by two fiber optic halogen lights 
that created PAR levels in the chamber similar to those in the 
mesocosms (~600 umol quanta m s ). Following photosyn- 
thesis measurements, samples were returned to their respective 
treatments and respiration estimates were made on the same 
samples at night between 1900 and 2300. During respiration 
measurements the chamber was kept dark using an opaque 
shroud. For each incubation (photosynthesis and respiration), the 
chamber was filled completely with seawater from the appropriate 
treatment and new seawater was used for each incubation. 
Seawater from the ambient treatment was used for both ambient 
and variable treatment incubations (because samples in the 
variable treatment were not placed in the elevated treatment until 
0000), and seawater from the elevated treatment was used for the 
elevated incubations. Oxygen concentrations were recorded each 
minute and photosynthesis and respiration rates were calculated 
from the linear slopes over time. Metabolic rates were normalized 
to surface area of the coralline algal samples, estimated by image 
analysis (ImageJ, NIH) of digital photographs of each sample. Net 
calcification rates for each sample were calculated by converting 
the buoyant weight gained during the experiment to the dry 
weight gained using the density of calcite (2.71 g cm ^ [78] and 
normalizing the rates to surface area. 

Statistical Analyses 

Normalized rates of net calcification, gross photosynthesis, 
respiration, and net photosynthesis were analyzed using a mixed 
model ANOVA, with habitat origin (upstream/ downstream) and 
pC0 2 treatment (ambient, elevated, variable pC0 2 ) as fixed 
factors and tank as a random nested factor. If the tank factor was 
not significant in the initial analysis (P>0.25), it was eliminated 
from the subsequent analysis and a two-way fixed factor ANOVA 
was used [79]. Additionally, an explicit a priori contrast tested 
specifically for the effect of the variable pC0 2 treatment on 
coralline algae from upstream and downstream habitats. The 
assumptions of normality and homoscedasticity were examined 
prior to all analyses through the examination of residuals and 
Cochran's test, respectively. 
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Figure 1. Mean daily and hourly pC0 2 of treatment conditions. 

Mean (± SE) daily pC0 2 of ambient (open circles) treatment tanks 
(N = 3), elevated (closed circles) treatment tanks (N = 3), and the overall 
treatment means for the duration of the experiment (dashed line) (A). 
Mean pC0 2 of ambient and elevated treatments measured over one 
diel cycle (B). The solid line shows the transition of the variable 
treatment samples from ambient pC0 2 into the elevated pC0 2 
treatment 0000 and then back into the ambient treatment at 0600. 
pC0 2 was calculated using C0 2 SYS and measurements of total alkalinity 
and pH (total scale). 
doi:1 0.1 371 /journal.pone.0087678.g001 

Results 

Response Variables 

The nested tank effect in the mesocosm experiment was not 
significant (P>0.25) in the initial analysis of net calcification so it 
was removed from the subsequent analysis [79]. Net calcification 



Table 2. Daily integrated pC0 2 
downstream backreef habitats. 


for upstream and 




Habitat 


Diel Period 


pC0 2 (uatm) 


Upstream 


Day 


424.8 ±9.4 




Night 


421.0±9.9 


Downstream 


Day 


385.8±5.9 




Night 


700.4±55.5 



Values (means ± SE) of in situ pC0 2 from daily integrated water samples at 
upstream and downstream sites of collection (N = 5 days). Upstream sites were 
located at the reef crest and downstream sites were located 200 m shoreward 
of the reef crest. Continuous water samples were collected simultaneously and 
autonomously from each site during the day (from 0800 to 1800) and during 
the night (from 1900 to 0400). 
doi:1 0.1 371 /journal.pone.0087678.t002 
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Figure 2. Calcification response of Porolithon onkodes to ocean 
acidification treatments. Mean (±SE) rates of net calcification for P. 
onkodes collected from upstream and downstream habitats after 14 
days of exposure to three pC0 2 treatments. The asterisk denotes a 
significant difference of an a priori comparison between upstream and 
downstream origins under the variable pC0 2 treatment. 
doi:10.1371/journal.pone.0087678.g002 

of P. onkodes was significantly affected by an interaction between 
the pCO z treatment and habitat origin (F 2j42 = 4.02, P = 0.025), 
with significant effects of both pC0 2 treatment (F 2 42 = 18.42, P< 
0.0001), and habitat origin (F 142 = 5.18, P = 0.028) (Figure 2, 
Table SI). As demonstrated in previous studies, net calcification 
was reduced by over 70% in the elevated pC0 2 treatment 
compared to ambient conditions and this trend was exhibited by 
individuals from both upstream and downstream habitats. The 
overall response of P. onkodes net calcification to the variable pC0 2 
treatment was an intermediate decrease that might be explained 
most parsimoniously by the proportion of time spent in the 
ambient (75%) and elevated (25%) pC0 2 treatments, although the 
overall decrease in calcification was closer to 35%. Notably, the a 
priori comparison between P. onkodes collected from upstream and 
downstream habitats revealed significantly different calcification 
responses to the variable pC0 2 treatment (F 1I4 = 22.43, P< 
0.001). Following exposure to the variable pC0 2 treatment 
corallines from the downstream habitats experienced a 17% 
decrease in net calcification as compared to the ambient controls, 
while corallines from the upstream habitats experienced a 41% 
decrease in calcification compared to ambient controls (Figure 2) 
There was no significant effect of nesting tank within treatment 
in the initial analyses of the metabolic response variables (P>0.25) 
therefore the tank effect was pooled in subsequent analyses [79] 
Rates of gross photosynthesis were not affected significantly by 
pC0 2 treatment (F 218 = 3.03, P = 0.086), habitat origin 
(Fi, is = 0.630, P = 0.443), or the interaction between factors 
(F 2 ^i8 = 0.738, P = 0.499) (Fig. 3A, Table SI). Similarly, there 
was no effect of pC0 2 treatment (F 2;18 = 0.414, P= 0.667), habitat 
origin (F ljl8 = 1.26, P = 0.276), or the interaction between factors 
(F 2 ,is = 0.020, P = 0.980) on P. onkodes rates of respiration (Fig. 3B, 
Table SI). However, the net photosynthetic rate of P. onkodes was 
significandy lower in the variable pC0 2 treatment compared to 
the ambient treatment (F 218 = 5.59, P = 0.013). There was no 
effect of habitat origin (F ljl8 = 0.041, P = 0.843) or the interaction 
between factors (F 2 i 8 = 0.882, P= 0.431) on net photosynthesis 
(Fig. 3C, Table SI). Data are archived at the Moorea Coral Reef 
LTER data repository (http://mcr.lternet.edu/data/). 
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Figure 3. Metabolic responses of Porolithon onkodes to ocean 
acidification treatments. Mean (±SE) rates of gross photosynthesis 
(A), respiration (B), and net photosynthesis (C) of P. onkodes collected 
from upstream and downstream habitats after 14 days of exposure to 
three pC0 2 treatments. There were no significant differences among 
treatments or origins for gross photosynthesis. Small case letters 
represent results of post hoc multiple comparisons for net photosyn- 
thesis where treatments sharing letters are not significantly different. 
doi:1 0.1 371 /journal.pone.0087678.g003 

Discussion 

Our results demonstrate that the response of a coralline alga to 
variability in pC0 2 depends upon the conditions in the original 
habitat. Porolithon onkodes collected from a downstream habitat, 
with prior exposure to high pC0 2 variability, calcified significandy 
more in the variable pC0 2 treatment than individuals from the 
upstream habitat, with history of exposure to lower and diurnally 
stable pCO z . These findings support our hypothesis that exposure 
to oscillating, high pC0 2 within the natural range of variability, 
had a modulating effect on coralline calcification compared to the 
static high pC0 2 treatment. We suggest that because the original 



habitat of P. onkodes had a significant impact on the coralline 
calcification response to variable pC0 2 , a history of exposure to 
dynamic carbonate conditions may facilitate plasticity in organ- 
ismal responses to ocean acidification. Our findings emphasize the 
importance of incorporating oscillating pC0 2 levels, where 
appropriate, into laboratory experiments because prior condition- 
ing to high pC0 2 may confound calcification responses to ocean 
acidification. 

Corallines exposed to the static treatment of high C0 2 calcified 
>70% less than those in the ambient treatment, regardless of the 
origin of collection. This finding corroborates several recent 
studies that have shown static, high pC0 2 conditions to have 
severe consequences for crustose coralline calcification and growth 
[40,55-60,80]. The reduction in calcified biomass following 
exposure to high pCO z may be a result of decreased calcification 
and/ or increased dissolution associated with a lower SI. Although 
the high pC0 2 treatments did not reach undersaturation with 
respect to aragonite or calcite saturation state (SI <1), calcification 
becomes increasingly difficult and potentially energetically more 
costly as it approaches the saturation horizon (£1=1) [3]. The 
difference in the reduction of net calcification for upstream and 
downstream organisms between the variable treatment and 
ambient pC0 2 treatments (41% and 17%, respectively) indicates 
that P. onkodes from downstream habitats were able to mitigate 
some of the negative effects of elevated pC0 2 on calcification, 
and/or the promoting effects of higher pC0 2 on calcium 
carbonate dissolution. The response of a marine calcifier to 
oscillating pC0 2 in a laboratory ocean acidification experiment 
has been demonstrated in few other studies, thus far. Dufault et al. 
(2012) showed that coral recruits exposed to oscillating pC0 2 
(between ~450 and ~850 uatm) experienced enhanced growth 
and survival compared to recruits exposed to both a static ambient 
and high pC0 2 treatment [36]. Together these findings indicate 
that variability in carbonate chemistry should be incorporated into 
future OA experiments in order to enhance the ecological 
relevance of laboratory manipulations. 

The increase in dissolved C0 2 associated with ocean acidifica- 
tion and C0 2 enrichment has the potential to fertilize photosyn- 
thesis by increasing the amount of substrate available to the 
photosynthetic enzyme Ribulose-l,5-bisphosphate carboxylase 
oxygenase (RuBISCO) [81,82]. High pC0 2 has been shown to 
stimulate photosynthesis in a variety of marine primary producers 
including seagrasses, phytoplankton, and some macroalgae 
[39,81,82]. However, species-specific responses may be influenced 
by the presence and activity of carbon concentrating mechanisms 
(CCMs) and whether a primary producer is carbon limited. We 
hypothesized that C0 2 enrichment would stimulate photosynthet- 
ic rates in P. onkodes by providing more substrate for carbon 
fixation. We found no effect of pC0 2 treatment on rates of P. 
onkodes gross photosynthesis or respiration, but a significant effect 
of the pC0 2 treatment on net photosynthesis, with highest rates in 
the ambient pC0 2 treatment, and lowest rates in the variable 
pC0 2 treatment. Because net photosynthesis was determined 
based on gross photosynthesis and respiration rates, our finding 
may be a statistical result that has little biological basis. The role of 
CCMs in coralline algae is not well understood, and differences in 
CCM activity may have contributed to the mixed results found 
here and in the recent literature. For example, Anthony et al. 
(2008) found that exposure to extremely high pC0 2 conditions 
(pH 7.6 and 7.8) decreased P. onkodes net productivity by 50% to 
100%, respectively [55]. Furthermore, Martin et al. (2013) found 
that high pC0 2 (700 |J,atm) decreased net photosynthetic efficien- 
cy of the Mediterranean coralline alga Lithophyllum cabiochae, and 
similar to our study found high pC0 2 to have no effect on 
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coralline respiration [83]. Respiration measurements were con- 
ducted between 1900 and 2300 using water from the appropriate 
treatments. For the variable pC0 2 treatments, samples were 
incubated with ambient water because samples were not placed 
into high pCO z conditions until midnight (0000). This method- 
ology did not account for the response of the variable treatment 
samples to high pC0 2 and, therefore, may underestimate the 
impact of the variable pCO z treatment on rates of coralline 
respiration. 

We suggest that the effect of habitat origin on the P. onkodes 
calcification response to variable pC0 2 provides evidence of in situ 
acclimatization to naturally occurring variable pC0 2 . Although 
the terms 'acclimation' and 'acclimatization' often are used 
interchangeably, here we follow the terminology of Edmunds 
and Gates (2008), wherein acclimatization refers to the physio- 
logical compensation of an organism to a regime of natural co- 
occurring environmental stimuli [84] . The relatively short distance 
between sites (200 m) and lack of physical barriers makes the 
possibility of genetic differentiation, thus adaptation, between 
upstream and downstream coralline populations unlikely. If local 
adaptation were the underlying mechanism to the significance of 
the origin by pC0 2 treatment interaction, we would expect both 
upstream and downstream individuals to exhibit the same 
response to variable pC0 2 . Form & Reibesell (2012) were among 
the first to document an example of acclimation of a marine 
calcifier to acidified conditions [42] . The cold-water coral, Lophelia 
perfuse, acclimated to acidified conditions after six months of 
exposure in a controlled laboratory mesocosm [42] . This provides 
an example of temporal acclimation to conditions similar to ocean 
acidification. In the current study, we suggest that corallines from 
downstream habitats are acclimatized to high pC0 2 , but only in 
the context of oscillating conditions that simulated the natural 
environment. 

It is difficult to discern the exact mechanisms driving lower 
coralline calcification rates in the static and variable pC0 2 
treatments, in part because calcification in coralline algae is not 
well understood. Corallines calcify intracellularly, where high-Mg 
calcite is deposited directly within the organic matrix of the cell 
wall [47,85]. In order for intracellular CaC0 3 precipitation to 
occur, the internal microenvironment must be supersaturated with 
respect to CaC0 3 (high pH). Calcification in corallines may be 
either biologically controlled or biologically induced [85] and 
influenced by both photosynthesis and calcification because they 
are inextricably linked. Biological control over calcification may be 
driven by photosynthesis, where photosynthesis directly provisions 
the energy necessary for calcification. Calcification may also be 
biologically induced as photosynthesis indirectly maintains the 
chemical conditions needed for CaC0 3 precipitation by producing 
OH~ ions [86]. For example, precipitation of magnesite and 
dolomite in P. onkodes cells is thought to be biologically induced by 
changes in pH due to metabolic activity [85]. Calcification is also 
influenced by external seawater chemistry, because calcified algae 
must uptake calcium ions (Ca + ) and carbonate ions (CO3 ) 
directly from the seawater or through the intracellular conversion 
of HC0 3 ~ to C0 3 2_ [87,88]. Biomineralization becomes more 
difficult when the CaC0 3 saturation state (Q) approaches under- 
saturation (< 1) [3], as is predicted to occur in the near future with 
ocean acidification. 

Plasticity in the coralline response to high pC0 2 variability 
therefore may have been driven by mechanisms associated with 
both photosynthesis and calcification. Identifying the mechanisms 
that promote coralline acclimatization to variable pC0 2 was 
beyond the scope of this study. Here we propose potential 
biological and mineralogical mechanisms that may have contrib- 



uted to the differences in coralline calcification by habitat. Two 
mechanisms (assuming biological control over calcification) by 
which corallines may maintain intracellular pH at the site of 
calcification, thus facilitating acclimatization to variable pC0 2 , 
are: 1) enhancement of photosynthesis from increased dissolved 
C0 2 concentrations, and 2) opportunistic use of HC0 3 — to 
mitigate the decrease in C0 3 2 . Photosynthesis may be fertilized 
by high C0 2 and thus facilitate supersaturation of the intracellular 
microenvironment. However, since photosynthesis was not higher 
in the variable pC0 2 treatment, this mechanism likely does not 
explain our calcification results. The second proposed mechanism, 
bicarbonate usage, applies to algae such as crustose corallines that 
can utilize HC0 3 ~ [88] through the conversion of HC0 3 ~ to 
C0 2 with the enzyme carbonic anhydrase. The enzyme-catalyzed 
conversion releases a hydroxyl ion (OH ) and increases the local 
pH environment in favor of CaC0 3 precipitation [54]. CCMs, 
including carbonic anhydrases (CA) and HC0 3 transporters, are 
important for supplying carbon to photosynthesis and calcification 
[89,90]. Production and expression of CCMs is energetically costiy 
and are influenced by environmental conditions [90,91]. The 
temperate coralline alga, Corallina officinalis, exhibited a trend of 
increasing internal CA activity following exposure to high pC0 2 
[92], and an inverse trend between external CA activity and C0 2 
concentration [93]. Corallines frequently exposed to high pC0 2 
therefore may regulate CCMs, such as CA activity, in order to 
mitigate decreased CO3 2 concentrations associated with OA by 
switching to HC0 3 ~ utilization [88]. A third mechanism 
(assuming abiotic control over calcification) may have been 
differences in calcification due to mineralogy. Coralline crusts 
from different habitats may exhibit differences in magnesium 
content, and thus corallines from upstream and downstream 
habitats may have differed in mineralogy from the start of the 
experiment. Mineralogy, rather than acclimatization, may have 
driven the calcification response to variable pC0 2 if abiotic 
calcification (cell infilling) by magnesite and dolomite [85] was 
affected differentially by pC0 2 conditions. In order to conclusively 
differentiate between mineralogical and physiological mechanisms, 
future studies should quantify coralline mineralogy before and 
after exposure to treatment conditions. 

Based on our findings, biogenic calcifiers may be found in 
greater abundances at microhabitats with lower daily pC0 2 
variability, a finding that corroborates field-based studies that 
showed the distribution and abundance of benthic calcifiers to be 
correlated with pC0 2 /pH variability [12-14,16]. As seen in other 
studies of shallow reef flats [7,12], our data show that shallow reef 
environments may be shaped by diel pH variability. Our findings 
are corroborated by seawater chemistry measurements from a 
shallow reef system in Okinawa, Japan, where the alkalinity and 
DIC are consistendy lower on the back reef during the day than on 
the reef crest [10]. Furthermore, the Tiahura reef transect, a 
similar section of back reef 6.5 km west of our sites on the north 
shore of Moorea, has been revisited since the early 1980's and 
calcification and carbonate chemistry have been thoroughly 
characterized, supporting our characterization of the upstream 
and downstream sites as having high and low pC0 2 variability 
respectively [31,67,94-97]. In addition to differences in the 
magnitude of carbonate chemistry variability, other factors that 
may influence coralline growth rates within a habitat include 
grazing intensity, nutrient availability, water flow, and light levels. 
Our coralline samples were collected from the same depth at 
upstream and downstream sites and samples were acclimated to 
the same light conditions prior to the start of the experiment. 
Previous light history in the original habitat (at the same depths) 
likely did not impact the observed growth rates among treatments. 
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Although upstream and downstream habitats may be character- 
ized by different flow regimes, this experiment was conducted 
during the Austral winter (May-June), when differences in flow 
between the upstream and downstream site are minimal [66]. 
Furthermore, flow-through water within mesocosm tanks, aera- 
tion, and water pumped from the chilling system created water 
flow within the tanks (~10 cm s -1 ). Although many factors 
influence coralline growth, the significant results from our 
laboratory experiment and the previous literature suggest that 
changes in the patterns of calcification rate among treatments were 
driven most likely by the carbonate chemistry conditions in the 
original habitats. 

This study documents a modulated response of P. onkodes 
calcification to high pC0 2 in the oscillating OA treatment which 
was based upon the natural range of pCO z experienced in situ. The 
significant impact of the original habitat on the coralline 
calcification response to pC0 2 treatments indicates that variability 
in carbonate chemistry likely is an important factor influencing 
biogenic calcification rates. These results are important because 
they emphasize that the original environmental conditions of 
individuals used in OA experiments can have a significant effect 
on experimental outcomes. Variability in carbonate chemistry in 
the original habitat may be a potential confounding factor in past 
OA studies, and may explain some of the ambiguity in calcification 
results in the literature. Further, we suggest that organisms in 
highly variable pC0 2 environments may be acclimatized to near- 
future changes associated with ocean acidification. Studies should 
continue to explore the adaptive potential of marine organisms to 
the changing environment by exploring both in situ acclimatization 
and adaptation. Our findings have important implications 
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